Systematics of heavy quarkonia from Regge trajectories on (n, M 2 ) and (M 2 , J) planes 
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In this paper we show that heavy quarckonium states, similar to light mesons, form Regge tra- 
jectories in (n, M 2 ) and (M 2 , J) planes and the slope of these trajectories is independent on the 
quantum numbers of the mesons. This fact can be useful for the prediction of the masses of heavy 
quarkonia and the determination of the quantum numbers of the newly discovered states. 



1. INTRODUCTION 

It is well known that for wide range of potentials the only singularities of the non-relativistic scattering amplitude 
are the poles in the the orbital momentum I plane (Regge-poles). The position of this poles depends on the energy, 
i.e. / = ct(s), where a(s) is the Regge trajectory that corresponds to physical particles and resonances. 

The collection of the values of the masses of light mesons qq is described fairly good by the linear Regge trajectories 
both on (n,M 2 ) and (M 2 , J) planes (here M is the mass of the meson and n and J are its radial quantum number 
and total spin respectively). The linearity of the dependence 

1 9 1 9 

n = —M 2 + n , J = — M + Jo 
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with the universal slope coefficient fi 2 ~ 1.2 GcV 2 is typical for the "string" nature of the mesons. In the series of 
works 0, the attempts were made to describe the relativistic string with fermions at the ends. In the limit of 
massless quarks the relativistic string model leads to linear Regge trajectories J = a'M 2 + a(0). When the quark 
masses are taken into account the trajectory deviates from linear for lower states. Another reason for such deviation 
is the screening of the quark-antiquark potential because of the additional gg-pair production. This effect results is 
important for higher states of the mesons. 

The spectrum of light mesons was analyzed in the works 0, Q and the Regge trajectories for this states were 
presented. Below we will perform the same analysis for heavy quarkonium states using the whole experimental 
information available at the moment. Such an analysis could be useful both for the comparison of the experimental 
data with potential model predictions and for the determination of the parameters of regge trajectories for heavy 
quarkonium. 

2. CHARMONIUM 

The most extensively studied charmonium mesons are vector states (J PC — l -- )- For this states the value of the 
total spin of the quark-antiquark pair is S = 1, while for the orbital momentum of this pair the values L — and L = 2 
are allowed, so the mixing of the states with different values of L is possible. This effect, however, will not change 
the masses of the particles significantly and we will neglect it in our paper. On fig^i,b we show the experimental 
values of the masses of vector charmonia with L = 0, 2 (red stars) and the predictions presented in the works Q and 
(they are labeled by the symbols ■ and A respectively). It is clearly seen that the experimental values are fairly 
well described by linear Regge trajectories and the mass of the meson M and its radial quantum number n satisfy 
the relation 

Ml = ,r,n ■ (1) 

Here the slope coefficient 

/4 = 3.2 GeV 2 (2) 
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is the same both for L = and L = 2, and the parameter Mq depends on the value of the orbital momentum. In the 
tabled we show the values of this parameter, the experimental values of vector meson masses and our predictions for 
the masses of excited mesons obtained with the help of the formula . 

The experimental information about the mesons with other values of L and/or S is more poor, since the observation 
of this states is more complicated. For each of the sets of the values of this quantum numbers the mass of at least 
one meson is, however, known, so we can determine the parameter Mq in the relation Q and using a known value 
of the slope coefficient (i 2 construct the corresponding Regge trajectory. In the table ^ we give the experimental 
quarkonia masses, the values of the parameter Mq and our predictions for the masses of the excited states. It should 
be mentioned the the slope coefficient /j? is universal for all trajectories and is defined according to the relation 
The trajectories are shown on the figures ^-g. 

On the figure^ 1 we show the Regge trajectories on the (J, M 2 ) plane for the ground states of charmonia 

J = a(M 2 ) = a'M 2 + a(0). (3) 
Here we used the value of the slope coefficient from the proceeding analysis: 



and the interceptions a(0) for different trajectories were found to be 



J/^(lS),Xc2(lP) 
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-2., 
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a(0) = 


-2.8 




a(0) = 


-3.5 



These values are almost in the range |a c (0)| > 3 4-4 that was obtained in the work with @ the help of QCD sum 
rules. The intercept of the parent trajectory was also presented in Q, where a(0) was expressed through the value of 
the i/)-meson at the origin. In that work the value a c (0) = —3.5 ± 0.6 was found and it agrees well with the results of 
the present paper. 

It is clearly seen that the positions of the physical states are described by the linear Regge trajectories with 
the universal slope with a pretty good accuracy. It should be mentioned that for the leading Regge trajectory the 
exchange degeneration holds, i.e. J/tjj and x C 2-mesons lie on one trajectory. This fact is not the result of the parameter 
fit, but automatic consequence of the correct description of the charmonia spectrum in (n, M 2 )-plane. It should be 
mentioned that the exchange degeneracy was used earlier for the determination of the intercept a(0), that is important 
for determination of c-quark wave function in J/ip, rj c , and D-mesons IqL Eo ). Ill f . 

3. NEW STATES 

Recently some new particles were observed. The production and decay channels of these particles indicate that they 
are excited charmonium states, so it would be interesting to consider the question of the disposition of this particles 
on the Regge trajectories presented in the previous section. 

3.1. X(3872) 

The first of such part icles was the X(3872)-meson. The properties of this state were widely discussed in the 
literature [H [H Q [H . 

X(3872) was observed by the Belle 0, CDFII [HI, D0 [3 and BaBar [HJ collaborations in the decay 

B ± -> J fC ± X(3872) -> K^tt'J/i/j, 

and the world average of its mass is 

M x = 3871.9 ± 0.5 MeV. 

The further experimental study showed that besides the X — * ir + it~ J/ip decay the decays X(3872) — ► jJ/ip |2Cj . 
X(3872) -> tt + tt-tt°J/%P H3 e X(3872) -> D°D 7r° [3 also exist and the branching fractions of these channels are 
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Table 1: The parameters of the Regge trajectories and masses of the mesons for (cc) sector. Our predictions for excited 
charmonium masses are shown in bold. 



linked by the relations 



Br^ -> K ± X)Br(X -> 7 V0 = (1.8 ± 0.6 ± 0.1) x 1(T 6 , 
Br^ -> K ± X)Br(X -> tt+tt^t/;) = (1.3 ± 0.3) x 10~ 5 , (4) 
Br(B ± -> K ± X)Bt(X -> D°D°n ) = (2.2 ± 0.7 ± 0.4) x 10~ 4 , 



Bi(X -> 7r+7r-7r» 
Br(X ^ tt+ti-VO 



1.0 ±0.4 ±0.3. 



For the total width of this state only the upper boundary is known: Fx < 2.3 MeV. If we assume that there are no 
other significant decay modes, then we can determine the branching fractions of the reactions listed above and set 
the upper bounds on their widths: 

Br « (7.4 ±0.4) x 10" 3 , T<(17±9)keV, 

Br « (5.3 ±0.8)%, T < (0.12 ± 0.02) MeV, (5) 
Br w (5.3 ±4.6)%, T < (0.12 ± 0.11) MeV, 
Br « (90 ±2.5)%, T < (2.1 ± 0.6) MeV. 
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Since the decay X(3872) — > 7^ is allowed, the charge parity of X(3872) should be positive. The angular distribution 
in the X(3872) — > 7r + 7r - ^ channel rules out the possibility of the scalar meson [2lj . 

Among the particles listed in the table n\xco{2P), Xci(2P) and T] C (3S) mesons have masses that are most close to 
the mass of X(3872). Since the case of the negative charge conjugation parity is forbidden by X(3872) — > yip decay 
and the case of the scalar meson contradicts the angular distributions, the only variant left is X(3872) = x c i(2P). 

There are, however, some arguments against this assignment. First of all, the upper bound for the width of the 
radiative decay T(X(3&72) — > 7^) < 17keV is less than theoretical predictions (for example, in [j| one can find the 
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values r[xci(2P) — > Jip] = 30 4- 60keV). Second argument is that the decay V(3872) — > pip — > 7r + 7r _ ^ implies that 
V(3872) is isovector, so it cannot be a charmonium. In the work |l3| an alternative model is considered. According to 
this work X(3872) is a loosely bound D°Z) *-molecule (deuson). The mass of V(3872) is surprisingly close to D a D°* 
threshold 

M D o + M D o, = 3871.2 ± 0.6 MeV 

and in the case of zero orbital momentum of the mesons in this molecule its quantum numbers should be equal to 
J PC = 1 ++ . This assumption explains pretty good both the isospin violation and the smallness of the radiative decay 
width. 

It should be mentioned that such an explanation have some serious drawbacks. First of all, the X(3872) mass 
is above D°D°* threshold. Secondly, the production probability of such a molecule should be smaller than the 
experimental value 

Bi(B + -> K + X) = (2.5 ± 1.0) x 10~ 4 . 

This value was obtained from the experimental results and JSJ and has the same order of magnitude as the 
production probabilities of other charmonium states in similar decays |l4j | . The smallness of the production rate is 
caused by the fact that the size of the loosely bound molecule should be lager than that of the (cc) state, and the 
value of its wave function at the origin should therefore be smaller j2^ . This problem could be avoided if we assume 
that X(3872) is a mixture of the charmonium and £>Z?*-molecule (in the weak decay of £?-meson (cc) component 
of this mixture is produced, and the isospin violation in the decay V(3872) — > pip is explained by the presence of 
Z).D*-component), but the other partner of V(3872) was never observed. 

We think that the other explanation looks more plausible. In the works 0,H3 it was shown that closeness of DD* 
thresholds (the mass of X a little bit larger than the D°D*° threshold and a little bit smaller than the D + D*~ one) 
increases the probability of the mixing of Xci(2P) with other charmonium states. Such a mixing leads to the decrease 
of the radiative decay width and can be a reason for the isospin violation. 

3.2. Z(3930), X(3940), Y(3940) 

The state Z(3930) was observed by the Belle collaboration in the reaction 77 — > DD [24[ and its mass equals 

M Z (s93o) = 3931 ±4±2MeV. 

The best candidates with the required charge parity from the particles listed in table are Xci(2P), r) c (3S) and 
X C 2(2P) mesons. The production and decay channels rule out first two variants, so we assign the quantum numbers 
jPC = 2 ++ tQ £(3900). 

The X(3940) particle was observed by the Belle collaboration recoiling against ip in e + e~ collision [2^. The mass 
of this particle is 

M X (3940) = 39 43 ± 6 ± 6 MeV 

and the dominant decay chanel is X(3940) — > D*D. The best variants for this state are Xci(2P), h c (2P) and rj c (3S) 
mesons. Production and decay channels exclude first two variants, so we choose the last assignment for X(3940). It 
is interesting to mention, that our prediction for rj c (3S) mass is only 19 MeV above the mass of V(3940), while the 
predictions of other models are 4040 — 4060, i.e. approximately 100 MeV too high [l3|. 

Belle has also observed one more particle in the region of 3940 MeV [2(| — V(3940). This meson was observed in 
the decay 

B -> if F(3940) -> KujiP 

and its mass is equal to 

M Y (3940) = 39 43 ± 11 ± 13 MeV. 

The decays V(3940) — > D^D have not been seen, so it is possible that V(3940) and V(3940) are dictinct states. 
Because of the existance of V(3940) — > uij} the charge parity of F(3940) should be positive. The best variants 
for this particle from table ^ are Xc/(2-P) mesons. Since the assigments for \ci 2(2P) are already chosen, we set 
V(3940) = Xco(2P). 
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4. BOTTOMONIUM 



In the previous sections we have shown that the charmonium states are described fairly good by the linear Regge 
trajectories on (n, M 2 ) and (M 2 , J) planes with the slope coefficient 

/4 = 3.2 GeV 2 , a' =4-- 

Mcc 

The similar analysis for bottomonium states (i.e. the (bb) mesons) shows that in this case the linear Regge trajectories 
describe the mass spectrum much worser. From figure [5^,, where all known vector bottomonia with the value of the 
orbital momentum L = are shown, one can see, that the slope of the Regge trajectory is not constant, but decreases 
with the increase of the radial quantum number n, so the linear Regge trajectory contradicts the experimental data. 
This fact should not be a big surprise. The linear Regge trajectories are typical for the "string" model of the quark- 
antiquark interaction under the assumption the the mass of the quarks tends to zero. Such an assumption is justified 
well for light mesons and agrees well with the experimental results. In the previous sections it was shown that the 
Regge trajectories of charmonium states are also linear, although the slope of these trajectories is different from that 
in the light meson sector. In the case of (bb) mesons, however, we see the significant deviations form the linearity. 
The possible reason of such deviation is that the massless quark approximations fails in this case. In the series of 
works (for example jll l2l 1271 l28j) the results obtained in the framework of the string model with massive quarks are 
presented and these results are in qualitative agreement with the experimental picture. 

The same result can be obtained in the framework of the potential models. Let us consider the widely used Cornell 
potential of the quark-antiquark interaction 29] : 

V(r) = —— + nr + c. (6) 
3 r 

This potential combines the main characteristics of the quark interaction known from QCD. At small values of the 
distance between quark and antiquark r <C r$ ~ y/ a s /n ~ 0.4 fm the leading term of this expression is the first one and 
we observe the asymptotical freedom. For large distances, on the contrary, the second term gives the main contribution 
and we observe the confinement. For the light mesons and charmonia the second case holds, so we can neglect the 
coloumbic term in the expression © and obtain the linear Regge trajectories on (n, M 2 ) and ( J, M 2 )-planes 

M 2 = 8k \2n + J + + const (7) 

that are in good agreement with the experimental data. For the lower bottomonium states, on the contrary, the first 
case takes place. After neglecting the second term in the expression |J5J we get 



M 2 = ^ — - | - ^ + const. (8) 



64a 2 m 2 1 
~9 (n + J + 1) 

Linking the expressions J7J and (JSJ) one can obtain the interpolation formula [3 



M 2 = 8K (2n + J + l)- ¥ ^ TW + Ml 



that joins both two limits. For large values of the quantum numbers the trajectory turns into linear with the slope 

^bb 



= 16k, 



and for small values of the quantum numbers the spectrum of the nonrelativistic quarkonium is restored. 

On the figures |2Ji-f we show the Regge trajectories obtained with the help of the formula © for all bottomonium 
states. It should be mentioned that the slope coefficient 

= 4.1 GeV 2 

is universal for all this trajectories, while the parameters b and Mq depend on the spin and parities of the particle. It 
is clearly seen that the experimental states lie good on this trajectories and deviate from linear ones for small values 
of the radial quantum number n. The values of these parameters, as well as known experimental masses masses and 
the predictions obtained with the help of the formula © are presented in the table [21 
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Table 2: Regge trajectories parameters and masses for (bb)-mesons 
. Our predictions for excited bottomonium masses are shown in bold. 



Up to now we have used the formula © for the construction of the Regge trajectories on (n, M 2 ) plane and the spin 
of the particle for each trajectory was fixed. This formula can be used also for the construction of the trajectories on 
(M 2 , J) plane (Chew-Frautchi plot). To do this one needs to fix the radial quantum number n and solve the equation 
@ for the spin of the particle. There are three solutions of this equation, but only one of them is physically sensible 
(the others give negative values for the real part of J and we will not consider them here). As it was shown earlier, 
in the case of charmonium states the parameters of the Regge trajectories on (n,M 2 ) and (M 2 , J) planes coincide, 
so we have checked this property for bottomonia. The trajectories for lightest bottomonium mesons (i.e. n = 1) 
are shown on figure [2h and one can see that they agree well with the experimental data. One can observe that the 
exchange degenerations holds for the 1 and 2 ++ states, as it was in (cc) case. We did not perform any fits to 
obtain the values for the parameters of these trajectories. Instead of this we used the parameters for S wave vector 
bottomonia (fig.^) for the upper curve, and 1++ and ++ bottomonia (figures [2Ji and [2:) for the middle and lower 
curves respectively. In the limit of large M 2 all these trajectories become linear and the interceptions form them are 
equal to 



T(1S), XM (1P) 
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-0.88 


% (LS),X6i(LP) 


a(Q) = 


-1.2, 
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These estimates, unfortunatelly, are not reliable. According to the QCD-sum rules restrictions 6] for the parent Regge 
trajectory we have 

a b (0) < -7 ±-8. 

The value of the same order of magnitude was presented in 0, where the intersept a(,(0) was expressed though the 




Figure 2: Regge trajectories for bottomonium states. — experimental values, A — prediction presented in tj ■ 
results of Q, D — our predictions 
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leptonic width of the T(lS')-meson: 

05(0) = -11.5 ±0.5. 

There is also another argument against our values of the interceptions. It seems logical that the absolute value of the 
intercept of the corresponding Regge trajectory should be larger for larger quark masses. 

In the work [3l| other parametrization for Regge trajectories was used, and at first sight this parametrization is 
free from the mentioned drawback: 

J = a(M 2 ) = a(0) + 7 (Vf - - AP \ , (10) 

where the "slope" 

7 = 3.65 ±0.05 GeV -1 

is universal and the values of T and a(0) are 

a(0) = -14.81 ±0.35, 
VT = 12.50 ± 0.02 GeV. 

It is clearly seen that the behaviour of this trajectory is completely different from the one considered in our paper. 
First of all, according to the expression the mass of the meson cannot exceed the threshold VT. The second 
difference is that the trajectory presented in [3l| is linear for small meson masses and deviates from hear near the 
threshold, wihle the trajectory used in our pape is, on the contrary, linear for large M 2 and deiviates from linear for 
small. 

The expression (llOf) was obtained in the framework of the generalized string model, where the tension of the string 
depends on the quark separation: 

a(r) 



1 + (ir/ir) 2 ' 

and the parameters T, <7 an d A* are connected by the expression T = (Xo/fj,. The same result can also be obtained in 
the framework of potential models if one use the potential 

a r, r < 1/fx, 



V B r(r 



^arctan^r)- ' (H) 

7T/X cr /(2ju), r > 1/fi. 

The asymptotical behaviour of this potential is completely different from the behaviour if the Cornell potential used in 
our work. For example for small values of quark separation Vb t (j') is linear with respect to r, wile for Cornell potential 
the coulomb term dominates. Thus, there is no assimptotic freedom in (|1 1|) . Other difference is that for large r the 
Cornell potential increases infinitely, wihle T^Br( r ) tendes to a finite constant. This could be tghe reason for the meson 
mass threshold existance in l31l. 



5. CONCLUSION 



The new states discovered recently in the (cc) sector open the question of their classification and the reliable 
prediction of heavy quarkonia masses. Such pred ictions were obtained in the framework of different potential or 
lattice models (see for example 0, 13^, H3, l34l l35l 136)1. but the results of this calculations depend strongly on the 
choice of the model parameters. For the mass of the Xco(2-P) meson, for example, one can find in the literature the 
values from 3.822 GeV to 4.080 GeV |33j. Since there is no reason to prefer some prediction it seems important to 
have some independent criterion that can help to choose the right value. 

We think that this criterion could be the positioning of the mesons to the respective Regge trajectories on (n, M 2 ) 
and (J, M 2 ) planes (here n is the radial quantum number of the meson and J and M are its spin and mass). It is well 
known that the masses of the lights mesons can be described with a pretty good accuracy by the linear trajectories. 
In this paper we show, that this is, with minor changes, valid also for heavy quarkonia. Namely, the charmonia Regge 
trajectories are the straight line with the slope the is common for all charmonium mesons. In this paper we have used 
this trajectories to position to them recently discovered particles V(3872), Z(3930), V(3940) and V(3940). 

In the (bb) sector the situation slightly changes. For the small values of bottomonia masses we observe a deviation 
of the trajectories from linear, while for excited states linearity restores. This is the behavior that one should expect 
from different theoretical speculations, for example form the string picture of quark-antiquark interaction or potential 
models with Cornell potential. The results presented in our paper could be useful to improve the accuracy of (bb) 
masses predictions. 
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